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Abstract 

In order to draw out the essential behavior of the universe, investigations of 
early universe cosmology often reduce the complex system to a simple inte- 
grable system. Inflationary models are of this kind as they focus on simple 
scalar field scenarios with correspondingly simple dynamics. However, we can 
be assured that the universe is crowded with many interacting fields of which 
the inflaton is but one. As we describe, the nonlinear nature of these inter- 
actions can result in a complex, chaotic evolution of the universe. Here we 
illustrate how chaotic effects can arise even in basic models such as homoge- 
neous, isotropic universes with two scalar fields. We find inflating universes 
which act as attractors in the space of initial conditions. These universes 
display chaotic transients in their early evolution. The chaotic character is 
reflected by the fractal border to the basin of attraction. The broader impli- 
cations are likely to be felt in the process of reheating as well as in the nature 
of the cosmic background radiation. 

05.45.+b, 95.10.E, 98.80.Cq, 98.80.Hw 
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I. INTRODUCTION 



The inflationary paradigm strives to deliver a smooth universe from random initial con- 
ditions. If inflation is a robust attractor in the space of initial conditions, then it earns 
its claim to naturalness and genericity A universe which hosts many different fields, 

including an inflaton candidate, can develop an extreme sensitivity to initial values. This 
sensitivity marks the onset of chaos. Chaotic dynamics does not in itself destroy the robust- 
ness of an inflationary phase. Rather, it can lead to some powerful and perhaps observable 
implications for a realistic universe. For instance, a fractal pattern in the spectrum of den- 
sity fluctuations could be generated. Also, the final phase of inflation marked by reheating 
would unavoidably be a setting for chaos. 

In simple cosmologies, the ultimate fate of the universe can be predicted once a set of 
initial conditions is prescribed. In a closed cosmology for instance, it can be determined 
from the initial prescription if the universe inflates or collapses. A plot in phase space will 
show regions or basins within which all of the initial conditions lead to the same outcome. 
There will be basins of inflation and basins of collapse. If the dynamics is not chaotic, these 
basins of attraction^ are distinctly separated by smooth, regular boundaries. If the dynamics 
becomes chaotic, then the smooth boundaries begin to break up, ultimately becoming fractal. 

The models described in the following sections are chosen on the grounds of simplicity, 
and do not necessarily conform to any standard inflationary scenario. In our current models 
the primary inflaton is weakly coupled and essentially not dynamical, leaving the chaotic 
dynamics to the other scalar fields. Multiply coupled fields are quite natural in any particle 
theory. The Higgs field for instance must couple to the standard model fermions in order to 
induce fermion masses. In supersymmetric theories, a glut of coupled particles is expected. 
The fields behave effectively like nonlinearly coupled harmonic oscillators and so naturally 
bring on chaos. In future studies we intend to connect more concretely to specific hybrid 
inflationary models || and draw out implications for the spectrum of fluctuations or the 
end of a realistic model. The main aim of this paper is to illustrate the complex dynamics 
that can arise in relatively simple cosmologies. 

In addition to highlighting the appearance of chaos in inflationary cosmology, we aim to 
demonstrate the power of fractals as a quantitative measure of chaos in relativity. In general 
relativity coordinate independent measures of chaos are of vital importance. One of the most 
valuable measures of a chaotic system in flat space, the Lyapunov exponents, can be removed 
by a simple coordinate transformation in curved space. Thus the usual coordinate dependent 
measures of chaos become ambiguous in a relativistic context. Topological signals such as 
fractal basins, cantori or stochastic layers in phase space are needed for conclusive evidence 
of chaos. In order to search for cantori or stochastic layers it is necessary to construct slices 
through phase space known as Poincare sections. Since Poincare sections rely on quasi- 
periodic behaviour, the system must complete many cycles for a useful picture to emerge. 



Throughout the paper we loosely refer to attractors in phase space. In more formal terminology, 
the word attractor is reserved for dissipative systems. Strictly speaking, our attractors are just 
asymptotic regions of phase space. 
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Oftentimes relativistic systems are not so obliging as the evolution may end at singularities, 
such as the big crunch or inside a black hole. In these cases the dynamics is better suited 
to an outcomes based approach such as the study of fractal basin boundaries [|J. 

In a chaotic system the different possible outcomes will each have a basin of attraction 
in the space of initial conditions, with the basins separated by a fractal border. A specific 
examination of phase space will require a coordinate system to be chosen. One might 
worry that in a different time slicing, the fractal would disappear. This is not possible. A 
coordinate transformation must be smooth and differentiable. No smooth map can undo a 
truly fractal pattern as fractals are non-differentiable. While the features of a fractal may 
be altered by a coordinate transformation, the existence of a the fractal is unambiguous. 

The cosmological context we employ allows us to demonstrate the technique of fractal 
basin boundaries. What emerges is a definitive manifestation of chaos in cosmology. We have 
already commented that we focus on multi-field systems. If only one scalar field is present 
and the universe is closed, then there is still the possibility of chaotic regions in phase 
space . However, in order to generate chaotic dynamics with just one matter field, the 
universe must oscillate between expansion and contraction many times. The requirement of 
many bounces makes these otherwise interesting solutions unlikely if not truly unphysical. 
We consider single field scenarios in §JTJ. Since these bouncing candidates do not represent 
viable cosmologies, we turn our attention to many field systems and the demonstration of 
fractal basin boundaries in §|TV|-|V|. Also note, that while we study closed cosmologies, the 
chaotic transients can be seen in a universe which never collapses or bounces. When fields 
interact, the chaotic nature is therefore not limited to closed cosmologies. It is thus possible 
that there was a transient chaotic epoch in the history of our own universe. 

II. THE COSMOLOGICAL MODEL 

In the coming examples we consider closed Friedman-Robertson- Walker (FRW) uni- 
verses. For potential driven inflation to be successful the inflationary potential needs to 
be fairly constant. For our purposes the inflaton can be modeled by a simple cosmological 
constant. If this were the complete system there would be of course no chaotic dynamics. 
However, the universe is created bursting with matter fields. We model the matter content 
by a variety of conformally and minimally coupled fields. These matter fields interact and 
can incite chaos. For the inflationary cases at hand then, the chaotic behaviour is principally 
matter driven, i.e. chaos in T^ v causing chaotic evolution of g^ v . 

The physical picture is that of a closed, pre-inflationary universe just exiting the Planck 
era. The space of initial conditions is probed by assigning three possible outcomes. Either 
the universe inflates forever, inflates for short spurts but then collapses, or collapses without 
any inflationary event. The outcome depends on the relative sizes of the various kinetic 
and potential energies in the matter fields. As the interaction between the matter fields is 
turned up the boundary which separates inflating from noninflating initial conditions blurs, 
eventually becoming fractal. 

While we only consider closed models in our outcomes based approach, the appearance 
of chaotic transients will be generic, regardless of curvature. This becomes clear since chaos 
is also nested within a given outcome basin. For instance, a universe will often go through 
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rocky beginnings, enduring many fits of inflation before taking off smoothly or collapsing. 
Within a collapse basin, the sensitivity to the initial conditions shows up as a random scatter 
in the maximum radius of the universe or in the final value of the fields. If the universe has 
managed to inflate by several e-folds there is no turning back as the kinetic energy which 
might interrupt inflation quickly redshifts away. This feature of the de Sitter attractor is 
often referred to as cosmic baldness. We shall see that while the de Sitter attractor might 
end up bald, it can have very hairy beginnings. 

We shall consider FRW universes described by the metric 

*' 2 - r 2 dsi 




where t is cosmic time and r is conformal time. Throughout, dots will denote derivatives 
with respect to cosmic time and dashes will denote derivatives with respect to conformal 
time. The matter Lagrangian will contain various combinations of conformally, and 
minimally, $, coupled scalar fields with a variety of interaction terms described by the 
potential V(^, $): 

c M = v - l -d^d^ - l -d^d^ + -^nv 2 . (2) 

For comparison with the conformally coupled term, i.e. the last term in eqn the gravi- 
tational Lagrangian is Cg = — fk^?.- We have chosen units where 4nG/3 — c — 1. In terms 
of cosmic time t the field equations read 

<£ + + d^V = , (3) 

^ + 3H^ + d^,V + ^ = , (4) 
6 



d + a (2$ 2 + + H^f + — * 2 + ^>d^V -2V\=Q. (5) 



The Ricci scalar can be related to the scale factor through 1Z/Q = a/a + H 2 + k/a 2 . The 
Hubble expansion factor is given by H = a /a. The total energy of the system is 



H = H 2 + ^-<b 2 -(4> + Hq) 2 -^ 2 -2V = 0. (6) 
a 2 a 2 



The constraint eqn (|6]) can be obtained directly from Einstein's field equations and represents 
the first integral of eqn (|). 

Since we are dealing with the entire universe, the system is necessarily conservative. We 
use eqn (|6]) to ensure that energy is in fact conserved. It is amusing to notice that if we 
isolate the matter sector, this subsystem looks dissipative. Energy is lost to the gravitational 
field through the friction terms ~ 3if <3>. We can see the effects of dissipation within the 
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larger context of the Hamiltonian system. For instance we can watch the matter trajectories 
shrink down onto an attractor as the volume in phase space is dissipated (cf §[V|). 

In some cases it is profitable to recast the field equations in terms of conformal time r, 
and the rescaled variables if) = a^, = a$ and U = a 4 V: 



a" 



" -—<f> + d <l> U = (7) 



a 



if)" + kif) + d^U = (8) 

a" + ka + -(<//- -<p) --UU- $dJJ) = . (9) 
a \ a J a 

The total energy of the system can be expressed as 

a 4 H = (a') 2 + ka 2 -((f)'- -A - ((if)') 2 + kip 2 ) - 2U = . (10) 



A universe is said to inflate if the scale factor a accelerates in terms of cosmic time, i.e. 
a > 0. The cosmic time acceleration is given by 



aa" - (a') 2 1 

a = 3 = ~3 

a 6 a 6 



i \ 2 
a x 



2U - ipd^U - (if)') 2 - hj) 1 - 2 U' <f) 



(11) 



Notice that conformally coupled fields tend not to contribute to a positive acceleration. In 
other words, even with a potential oc \& n , conformally coupled fields do not drive inflation 
(unless n < 2). The generalization of the above equations to describe two or more scalar 
fields of either type is direct. 

The field equations are invariant under the combined rescaling 

a — ► xa, t — > xt, V — > — , (12) 

x l 

where x is a constant. This freedom is removed when we set our length scale by choosing 
dimensionfull values for quantities such as masses and coupling constants. 

As we describe below, the asymptotic solutions are of two kinds. The universe eventually 
approaches a smooth de Sitter phase or it ultimately collapses into a big crunch. In highly 
simplified models the division between these two outcomes can be expressed as a simple 
partition in the space of initial conditions. However, we shall see that even in simple models 
with two interacting fields the division is no longer clean, as the boundaries separating the 
different outcomes are no longer smooth curves but fractals. 



A. From integrability to chaos 

When the various scalar fields are massless and do not interact the equations of motion 
can be integrated exactly and there is no chaos. Before launching into the chaotic dynamics, 
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we take a look in this subsection at the two asymptotic possibilities, the big crunch and 
de Sitter expansion, which will be the basis of our outcomes approach in the following 
sections. We also show the phase space portraits for the non-interacting, closed system. 
In Ref. 0, a detailed analysis was given of the phase space portraits for a single, massive, 
minimally coupled field in a universe with arbitrary curvature. According to their portraits, 
the trajectories drawn are untangled and therefore are not chaotic. As the authors noted, 
there do exist a set of measure zero oscillatory solutions which do show chaotic behavior . 



At the close of § pl| we return to discuss this special set of perpetually bouncing solutions. 

For now we demonstrate the non-chaotic phase space for a universe full of noninteracting, 
garden variety fields (both minimally and conformally coupled). Since it is impossible to tell 
one scalar field from another in the absence of interactions, the general case of N minimally 
coupled and M conformally coupled scalar fields reduces to a universe with one field of 
each type. For a universe with a scalar field of each type and cosmological constant A the 
equations simplify to 

ip = A cos r , (13) 



B 

a? 



(14) 



a =±yA 2 + 2a 4 A + ^-a 2 . (15) 

The phase space is divided by a separatrix into two classes of trajectories, those that termi- 
nate at the big crunch and those that inflate. The form of the solutions can easily be found 
in the neighborhood of these two geometrically distinct attractors. When the dynamics is 
dominated by potential terms, such as a cosmological constant A, the universe undergoes 
exponential expansion and matter fields are redshifted away: 

a r-j e^ 2 ^ 1 ~ — ( — — J , 4> ~ a~ 2 , ip ~ cos r . (16) 



2A r, 

Conversely, when the dynamics is dominated by the kinetic energy of the matter fields or 
spatial curvature, the universe collapses to the big crunch at time r c : 

a ~ \t — Tel 1 / 2 , <p ~ a In a , ip ~ cos r (17) 

for B ^ and for B = 0: 

a ~ |r — t c \ , ~ constant , if) ~ cos r . (18) 

The separatrix that partitions these possibilities is defined by the trajectory with 

6{A 2 + VA* + 3B 2 ) 2 V ; 

For these simple, integrable cosmologies the basins of attraction for the big crunch and de 
Sitter attractors are separated by a smooth curve. This smooth curve is a portion of the 
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separatrix. In Fig. 1 we display phase space portraits in the (a, a') and (■0, ip') planes for a 
universe with B = and A = 1/8. The cross-hatched region is the basin of the big crunch 
attractor and the solid line is the separatrix. 

When interactions are included the separatrix breaks up and is replaced by a fractal curve. 
The gaps in the broken separatrix have the structure of a Cantor set. The broken separatrix 
no longer partitions phase space and trajectories may diffuse through it. For example, a 
universe that was destined to collapse in the integrable case might diffuse through the broken 
barrier and inflate. The breaking of the separatrix is reflected in the fractal nature of the 
basin boundaries for chaotic universes. The smooth basin boundaries shown in Fig. 1 should 
be compared to the fractal boundaries seen in Figs. 4 and 7. The break up of the separatrix 



is further described in 5 IV 



Even when interactions are included, the asymptotic behaviour of trajectories on either 
attractor is completely regular and non-chaotic. Examples of this fact are given in §[V]. In 
the parlance of dynamical systems theory, the attractors are neither strange nor chaotic. 
The chaotic behaviour is transient ||, and occurs when trajectories approach the broken 
separatrix. Physically this corresponds to an epoch in which the universe coasts with a 
fairly constant radius but with wildly varying acceleration. During this epoch the kinetic 
and potential energies in the system fight for supremacy and the universe teeters between 
collapse and violent expansion. If the kinetic energy wins the day, the universe collapses 
and the asymptotic solution can be found by neglecting all potential terms in the equations 
of motion. Conversely, if the potential energy wins the day the asymptotic solution can be 
found by neglecting kinetic energy terms. 

The transient nature of the chaos is similar to that found in the Mixmaster universe P,|10 



where it has been shown that the underlying attractors are neither strange nor chaotic |l I 
We remark that transient chaos appears to be the hallmark of relativistic systems. 



III. COSMOLOGIES WITH A SINGLE SCALAR FIELD 

Since we work with a closed FRW cosmology, there is only one parameter describing 
the gravitational sector - namely, the scale factor. The necessary elements for chaos are 
present if the scale factor interacts even with just one matter field. However, the dynamical 
timescale for the onset of chaos is longer than the life of one universe. The chaotic dynamics 
results as the two oscillators interact. Typically, at least a few oscillations are needed for 
the effects to surface. We discuss such an example in this section. On the other hand, if 
there are many interacting matter fields in the universe, then their chaotic evolution will 
make an impact during the lifetime of one universe. The examples of the following section 
reveal chaos on such short time scales. 

We begin with an example that is chaotic, but only on a timescale longer than the life 
of one universe. The model describes a single, conformally coupled scalar field in a closed 
(k = 1) universe. We choose the potential to have both a mass term and a cosmological 
constant A: 

U = ~m 2 aV + a 4 A . (20) 
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This example has previously been considered by Calzetta and El Hasi . The Hamiltonian 
takes the form 

cfH = (a') 2 + a 2 - ((^') 2 + ^ 2 + m 2 a 2 ij 2 ) - 2a 4 A = , (21) 

which, aside from the wrong sign for the gravitational contributions, is the Hamiltonian for 
two coupled harmonic oscillators. 

For a metric with only one dynamical degree of freedom it is always possible to perform 
a combined field redefinition and conformal transformation to a coordinate system in which 
the dynamics appears to be nonsingular. By using conformal, rather than cosmic time to 
describe the evolution of this system, the dynamical equations can be smoothly integrated 
past the big bang and big crunch singularities at a = 0. This allows many cosmic cycles to 
be considered if we continue the scale factor into negative values. When evolved through a 
series of cosmic cycles the system is clearly chaotic H|/]], as we might expect for nonlinearly 
coupled oscillators. It should be noted that the cosmic cycles are physically meaningless as 
all memory of the previous cycle is erased at each big crunch singularity. 

By introducing the fiction of cosmic cycles, the dynamics can be surveyed using the 
standard tools of Poincare sections (return maps) and Lyapunov exponents. Lyapunov ex- 
ponents measure the rate of separation of trajectories in phase space. Only if trajectories 
separate exponentially fast do they have positive exponents. Systems with positive Lya- 
punov exponents are said to exhibit sensitive dependence on initial conditions - one of the 
two ingredients of chaos (the other being the mixing and folding of trajectories). The in- 
verse of the positive Lyapunov exponents is referred to as the Lyapunov timescale. This 
timescale sets the dynamical timescale over which chaotic effects make themselves felt. In 
general relativity, Lyapunov exponents must be used with extreme care, if at all, as they are 
coordinate dependent. Indeed, a simple coordinate transformation can give a non-chaotic 
system positive exponents and a chaotic system vanishing exponents. 

Putting these reservations aside, we may compare the Lyapunov time to the time taken 
to complete a cosmic cycle, and infer whether or not chaotic effects can make themselves 
felt in the lifetime of a single universe. Typically, the Lyapunov timescale was found to 
be in the range 10 — > 1000 cosmic cycles. Even when the mass is taken to be very large, 
the Lyapunov timescale is always found to be greater than half a cosmic cycle, or in other 
words, the timescale for chaos to become important always exceeds the life of one universe. 
This result is easily understood. The chaotic behaviour is due to resonances between the 
two oscillating fields a and ip. In order for the resonance to take effect, both fields typically 
need to oscillate several times. However, a can only complete half an oscillation before the 
big crunch, making it exceedingly difficult for a chaotic resonance to occur. 

In Ref. || it was argued that chaos had been viewed within the span of one life cycle. 
Their conclusion was based on what appeared to be a scatter between initial values of the 
matter fields and the final values. The correlation between initial and final values of the 
scalar field ip was found to be 0.01. However, this low value for the correlation actually 
stems from a coarse sampling of a high frequency function. By regenerating Fig. 6. of Ref. 
H with a sampling rate that is ten times higher we see from Fig. 2. that the true correlation 
coefficient is 1.00. This confirms that the system shows no meaningful chaotic effects in the 
life of one universe. 
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Similar conclusions hold for universes inhabited by a single minimally coupled scalar 
field. Again, the equations of motion can lead to chaotic behaviour as they are nonlinear 
and have phase space dimension greater than 2. However, meaningful chaotic effects can 
only occur if the universe itself oscillates. The dynamics of an inflationary model driven by 
a minimally coupled scalar field with potential 

V = l m ^2 + ^^2_^2y ? (22) 



was studied by Belinskii et al. and with A = by Hawking |12| and Page ||. Typical 
trajectories were not chaotic. Rather, they see the universe smoothly evolve from the big 
bang to the big crunch with various amounts of inflation [|TJ. However, the inflationary 
potential allows some atypical trajectories for which the universe undergoes a number of 
nonsingular bounces [Q. Page [Q suggested that there exist an uncountably infinite but 
discrete set of perpetually bouncing universes with vanishing Lebesgue measure but non- 
vanishing fractal dimension. If Page's suggestion is correct, it would prove that the dynamics 
is chaotic as his "fractal set of perpetually bouncing universes" corresponds to what is now 
known as a strange repeller |TJ| . In contrast to the fictional cosmic cycles used to describe a 
conformally coupled scalar field, Page's bouncing universes are true, nonsingular solutions. 
However, these solutions have obvious drawbacks as plausible cosmologies. As remarked in 
Ref. HJ, the fine tuning required to arrive at these chaotic trajectories rules them out as a 
robust physical model displaying chaotic behaviour. Perhaps in a model of the early universe 
that generically displays nonsingular bounces we can hope to see interesting chaotic effects 
caused by an oscillating scale factor. In the absence of such a model we have to look to 
additional matter fields to provide the nonlinear resonances needed to incite chaos. 



IV. COSMOLOGIES WITH TWO CONFORMALLY COUPLED FIELDS 

If additional fields occupy the universe, then the scale factor will not be the principle 
source of chaos. Two scalar fields can oscillate many times in the lifetime of one uni- 
verse, leading to truly chaotic behaviour. To demonstrate the chaos we show the fractal 
basin boundaries for a universe which contains two conformally coupled fields which inter- 
act through the potential 

U = ^m\a 2 ^l + ^m 2 2 a 2 ^ 2 2 + X 2 ^ 2 + a 4 A . (23) 

The period of oscillation for each field is governed by its effective mass. We define the 
reduced effective mass for each field as the derivative with respect to the field of the field 
eqn (|^). In other words, M 2 has the form of d 2 W/dip 2 were W is anything which acts as a 
potential in the equations of motion: 

M a = (1 - m\^\ - rn^l - AAa 2 f 2 , (24) 
Mi = (I + mla 2 + 2\ 2 iPl) 1/2 , (25) 
M 2 = (1 + m\a 2 + 2\ 2 i) 2 ) 1/2 . (26) 
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Increasing mi, m 2 , A and A slows the recollapse of a and speeds the oscillation of ipi and fa, 
thus increasing the probability of chaotic resonances. However, if m\ or m 2 greatly exceed 
A, the resonances will be washed out and no chaos will be seen. Conversely, if m\ or m 2 
are both zero, the oscillations tend to freeze when ipi and fa hit small values, again making 
chaotic resonances unlikely. 

To gain some intuition we can find a simple analytic approximation which corresponds to 
a familiar chaotic system. During the majority of the universe's evolution, the scale factor 
varies much more slowly than the scalar fields so that a' /a <C fai/fa- The scalar fields behave 
like coupled nonlinear oscillators, adiabatically pumped by the slowly varying scale factor. 
To leading order we can ignore the adiabatic pumping all together and study the scalar 
field dynamics in a fixed background (a' rs 0). This approximation is particularly good 
for describing universes that are vacillating between collapse and inflationary expansion. 
Importantly, this is just the region where the chaotic transients occur that destroy the 
smooth separatrix of the integrable model described in § 11 A , When a' ~ the dynamics 
simplifies to that of two coupled oscillators: 

fa[ + + 2A 2 V^i = , (27) 

fai + u 2 2 fa + 2A 2 V^ 2 = , (28) 

where uj\ = 1 + mfa 2 is the fixed frequency of the uncoupled (A = 0) oscillators. The above 
system of equations describes a known chaotic system []i~4| , and the transition to chaos as 



A is increased can be studied using the Chirikov resonance overlap condition Having 
established that the fast variables fa an d fa behave chaotically, we can then consider how 
they backreact on the slow variable a. When looked at on timescales long compared to the 
periods of the scalar fields, the evolution of the scale factor is similar to Brownian motion, 
and can be described in terms of chaotic diffusion equations M. It is this buffeting of the 
scale factor by the matter fields that breaks the separatrix in the (a, a') plane and causes 
the universe to evolve in a chaotic manner. 

Returning to the full, unapproximated equations we numerically investigate the phase 
space of initial conditions. For a given set of initial conditions we can identify three main 
outcomes. The first possibility sees the universe expand and collapse without any inflationary 
burst. The second possibility sees the universe undergo one or many short bursts of inflation, 
but failing to become a macroscopic universe. The third possibility sees the universe sustain 
a prolonged and violent period of inflation resulting in the formation of a macroscopic 
universe. The first and second possibilities (coloured Black and Grey respectively) combine 
to form the big crunch basin of attraction. This artificial division of the big crunch basin is 
mostly for visual effect. There is a fourth possible outcome that should be mentioned. There 
are a set of trajectories with zero Lebesgue measure that oscillate eternally, never entirely 
collapsing or reaching the de Sitter attractor. These trajectories form the border between 
the big crunch and de Sitter basins of attraction. We will see that these trajectories belong 
to a fractal set of perpetually bouncing universes. In the parlance of dynamical systems 
theory, this set forms the stable manifold of a strange repeller ||. 

The three possibilities are displayed graphically in Fig. 3 for the choice of parameters 
(A = 0.0001, mi = 0,m 2 = 0.05, A = 1) and initial conditions {o(0) = 2, fa = 0.4, fa = 
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6, ip'2 = 20}. The initial values of ip[ are {—23.31, —23.32, —23.33}, and a'(0) is fixed by 
the Hamiltonian constraint. 

The fact that minute changes in the initial conditions can lead to such dramatic changes 
in the outcome suggest that the fate of our model universe is indeed chaotic. This suspicion 
can be confirmed by studying the boundary between the basins of attraction of the three 
outcomes. Since the basins are embedded in a six dimensional phase space, we are forced 
to consider lower dimensional slices through the boundary. In Fig. 4 we display a two- 
dimensional slice in the ipi - ip[ plane for universes with parameters and conditions identical 
to those used in Fig. 3. The three basins of attraction (Black, Grey, White) are dramatically 
intermixed strange basins, as at least a portion of the boundaries are fractal. The boundaries 
near the origin are regular and smooth while the outer boundaries appear fragmented. A 
detail of the outer region is shown in Fig. 5, visually confirming the fractal nature of the 
boundary. Repeated magnification reveals similar striated pictures on all scales. 

Rather than rely on these qualitative features, we may quantify the fractal nature of the 
boundary in terms of the fractal dimension. There are many definitions of fractal dimension 
that we may choose from, but the one best suited to our situation is the box counting 
dimension. On a two-dimensional slice through phase space we cover the fractal with a grid 
of squares with side length e. We then count the number, N(e), of squares needed to cover 
the fractal, i. e. the number of squares containing more than one colour. The box dimension 
ds is defined by 

, IniV 

d B = -hm- . 29 

For self-similar structures the formal limit e — > need not be taken, and in all practical 
situations we are only interested in the existence of such scaling laws over a large, but 
not necessarily infinite, range of scales. Since the fractal dimension is not invariant under 
homeomorphisms, it is not a true topological invariant. However, it is invariant under 
diffeomorphisms, so it does provide a topological measure in general relativity. The existence 
of fractal structures in phase space provides a coordinate independent signal of chaos in 
relativity. 

The importance of the fractal dimension of the basin boundaries can be described in 
terms of final state sensitivity [0 . Consider an initial configuration near the basin boundary, 
where the uncertainty in the initial conditions describes an N-dimensional ball of radius 5 in 
the N-dimensional phase space. The final state sensitivity f$ is the fraction of phase space 
volume which has an uncertain outcome due to the uncertainty in the initial conditions, and 
is given by 

f s = 6 a , a = N-d B . (30) 

For a non-chaotic system a = 1 and the final state sensitivity is directly proportional to 
the initial uncertainty. For chaotic systems however, < a < 1, and the uncertainty in the 
outcome is greater than the uncertainty in the initial conditions. For example, if a = 0.47, a 
50% reduction in the initial uncertainty only reduces the final state uncertainty by 28%. In 
this way, the dimension of the basin boundary is a direct measure of "sensitive dependence 
on initial conditions". 
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In Fig. 6 we display the plot used to determine the fractal dimension of Fig. 5. Because 
the three boundaries are densely interwoven in this case, we chose only to calculate the 
dimension of the boundary between the big crunch and de Sitter attractors, i.e. counting 
Grey and Black as one basin. Using an 840 x 840 grid we found the dimension to be 
1.58 ± 0.02. The grid size of 840 = 2 3 x3x5x7 was chosen as it has the most factors 
of any number below 1000. The curvature of the data points at small and large e is to be 
expected. For large e the covering is very inefficient, while for small e the squares saturate 
the resolution used to generate the fractal. These effects cause d B to tilt toward 2 for 
large e and toward 1 for small e. Despite these limitations, accurate fractal dimensions can 
be obtained very quickly and easily. For different choices of parameters we found fractal 
dimensions ranging from 1 to 1.96, essentially filling the allowed range d B = [1,2]. 

The boundary was found to be fractal on all possible two-dimensional slices. For example, 
in Fig. 7, the boundary is shown in the a - a' plane for a slice which intersects Fig. 3 along 
the line ipi = 1.0. The fractal dimension of this slice was found to be d B = 1-37 ± 0.02. 

While the previous chaotic pictures were typical of those found, the dynamics of the 
system is not always chaotic. For small values of A, mi, m-i (at fixed scaling x) the dynamics 
is near integrable and the basins are not strange, but regular. In Fig. 8 we increment A 
while keeping all other parameters and initial conditions fixed. The mixing of the basins 
is reminiscent of the blending of viscous fluids. The dimension of the basin boundary for 
A = 0.5 was found to be d B = 0.99 ± 0.02, which is consistent with a dimension of 1. 
So, within errors, this boundary is smooth and non-chaotic. To compare, the dimensions 
of the boundaries for A = 2.0 was d B = 1.16 ± 0.05 (Grey- White) and d B = 1.26 ± 0.05 
(Grey-Black). 

An important property of dynamical systems with strange attractor basins is that the 
chaotic dynamics is not restricted to phase space trajectories near the fractal boundaries. 
One way to see this might be to use the fiction of cosmic cycles to follow the evolution of 
trajectories starting in the big crunch basin. The Lyapunov exponents and Poincare sections 
for these trajectories would reveal chaotic behaviour across the basin. However, we are not 
really interested in effects which take longer than one universe's lifetime to make themselves 
felt. Instead we plot in Fig. 9 the correlation between the initial value of ipu and the value 
at the point of maximum expansion, ipi m . The graphs are for a ^ = slice through the 
big crunch basin of Fig. 8 with A = 2.0. The big crunch basin stretches from ipi — to 
ipi ~ 2.9 (and similarly for negative ipi). A general increase in frequency with increasing 
ipu requires that we use several plots, each covering half the region of the last, to cover the 
basin. Unlike the regular plot seen in Fig. 2, the relationship between initial and final values 
of ipi is highly erratic, with apparently random changes in frequency and amplitude. 

V. COSMOLOGIES WITH MINIMALLY AND CONFORMALLY COUPLED 

FIELDS 

The chaotic behaviour seen in the previous system is not restricted to conformally coupled 
fields. Similar behaviour is found for minimally coupled fields with the same choice of 
potential. The main difference in this case comes from the scalar fields themselves being 
a source of inflation, in addition to the cosmological constant. The acceleration in this 
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example is given by 



a = 2a (A + \ m \§\ + \m\<b\ + A 2 $ 2 $ 2 -¥ x - $ 2 ) , (31) 

where we have reverted to the unsealed field variables. These models are able to successfully 
inflate even when there is no cosmological constant, in a manner similar to Linde's "chaotic 
inflation" ||17|| . However, we did not see any strange basins when A = as successful 
inflation generally required the fields to become stuck high up in their potentials after just 
a few oscillations. Otherwise, their ability to climb high enough was lost due to friction and 
redshifting of kinetic energy. It may be that chaotic behaviour does occur when A = 0, but 
it is difficult to search for as the inflationary bursts must be followed for ~ 60 e-folds in 
comparison to the ~ 5 — 10 e-folds required to ensure we have reached the de Sitter attractor 
when A / 0. 

In Fig. 10 we display the basins of attraction in the $i - $i plane for universes with 
(mi = 0, m 2 = 0.04, A = 0.00005, A = 1.0) and fixed initial conditions {a(0) = 10.0, $ 2 = 
0.4, $2 = 0.16}. A detail of the outer boundary is shown in Fig. 11, where the dimension was 
found to be 1.54 ± 0.02. Because the scalar fields themselves contribute to the inflationary 
bursts there is typically far more of the Grey basin than we saw for conformally coupled 
fields. 

We can compare an analysis of the chaotic trajectories with the non-chaotic trajectories 
of §11 A| . While the universe is expanding (H > 0), we see from eqn (§) that the scalar field 
dynamics is effectively that of a damped harmonic oscillator. Conversely, as the universe 
contracts the dynamics is that of a pumped harmonic oscillator. This behaviour is apparent 
in Fig. 12, where we have displayed typical trajectories leading to the de Sitter and big 
crunch attractors. For the de Sitter attractor the scalar fields spiral into a fixed point as 
cosmic baldness asserts itself, while for the big crunch attractor the scalar fields first spiral in 
and then spiral out again as the universe collapses. The de Sitter attractor has an interesting 
structure when two fields are present as one scalar field gets locked at a constant value. The 
attractor is of the form 



a c exp(V2At) , (32) 
/ 3V2A \ 

$ 2 = *2 C exp —t cos ut , (33) 



(34) 




, / — / t 2u cos 2 ut + 3V2A cosa;t sintut" 
1 + A 2 $^ c exp(-3V2At) —j= + 



3^ 2cj(9A + 2u 2 ) 



u 2 = m 2 2 + 2A 2 $ 2 C - ^A . (35) 

For reference, the trajectory in Fig. 12a has $i c = 0.00573, $2c = 105.9 and to = 0.03796. 
The big crunch attractor is unchanged from the one-field case, and takes the form 
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a = a c (t c - t) 1/3 , 



(36) 



$1 ~ $ lc ln(t c -t), 



(37) 



$ 2 ~ $ 2c ]n(t c -t), 



(38) 



with $ 2 C + $| c = 1/9. The trajectory shown in Fig. 12b has t c = 200.004, a c = 13.68, 
$i c = 0.294 and $ 2c = 0.167. Since we are able to write down analytic solutions for 
trajectories on the attractors, it is clear that the de Sitter and big crunch attractors are 
non-chaotic. The chaotic behaviour seen during the evolution of the universe is restricted 
to the 5 — > 10 transient orbits seen in Fig. 12. The fractal nature of the attractor basin 
boundaries is due entirely to these brief chaotic transients. 

We close with a word on a mixed cosmology which contains one minimally coupled and 
one conformally coupled scalar field. The conformally coupled scalar field \1/ is taken to be 
massless. When ^'s coupling to the minimally coupled scalar field $ is small it behaves like 
radiation. The interaction potential is taken to be 



As before, the minimally coupled scalar field is able to contribute a negative pressure to that 
of the inflaton, thereby increasing the likelihood of inflationary bursts. Again, this increases 
the proportion of Grey over what we saw for two conformally coupled fields. 

By viewing the basins in the a - a plane we see some rather striking ink-blot and crystal 
boundaries. An example of this is shown in Fig. 13 for the choice of parameters (m = 
0.05, A = 0.0001, A = 2), and fixed initial conditions {$ = 0.2, * = 0.1, ^ = 0.08}. A detail 
of the Grey- White crystal boundary is shown in Fig. 14. The high degree of self-similarity 
of this fractal allowed a particularly accurate determination of the fractal dimension using a 
standard 840 x 840 grid. The dimension was found to be d B = 1.484 ± 0.005. The visually 
less fractal Grey-Black ink-blot boundary was found to have a smaller fractal dimension of 
d>B = 1-11 ± 0.05, with the larger error due to a lower degree of self-similarity. 

In addition to studying different combinations of minimally and conformally coupled 
scalar fields, we also considered a variety of polynomial potentials V(&, \1/). The qualitative 
results were the same for all cases, showing that chaotic evolution was a generic feature of 
all multi-field models. 



The early universe is likely to be home to many interacting fields. We have show that if 
these interactions are sufficiently strong, the evolution of the universe will be chaotic. The 
fractal basin boundaries reveal the chaos in a coordinate independent manner. Additionally, 
the method does not require one universe to pass through many cycles. We now have to 
ask how prevalent chaotic behaviour will be in particular models of inflation, and what 
implications it might have for processes such as reheating or galaxy formation. 



V = -m 2 $ 2 + \ 2 q 2 $ 2 + A. 



(39) 



VI. DISCUSSION 
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One model of inflation where chaotic dynamics is bound to be important is hybrid infla- 
tion. Hybrid models employ several interacting scalar fields, and arise naturally in various 
super symmetric theories where the breaking of large gauge groups employs many Higgs par- 



ticles [flq| . If in the early stages the universe inflates in jolts, an exciting possibility exists for 
the spectrum of primordial density fluctuations. Chaotic resonances could lead to a fractal 
power distribution, perhaps helping to explain the hierarchical clustering seen in the current 
universe. Moreover, chaotic evolution of the scale factor would leave a unique imprint on 
the gravitational waves produced during inflation ||19|| . Any chaotic behaviour would have 
to occur within the last ~ 60 e-folds of inflation to be observable today. This would require 
some artificial fine tuning in a single field model but may be more natural in a hybrid model. 

Even in inflationary models where chaotic evolution is unimportant at early stages, chaos 
is sure to play an important role at the end stages. To illustrate, consider again a hybrid 
model. The fields to which the inflaton couples dictate the occurrence of the true vacuum and 
so control the nature of the exit from inflation. The setting is prime for chaotic interactions 
which would certainly impact on the exit style. More generically, at the end of any inflation 
model, the universe reheats as the inflaton oscillates about the minimum of its potential. 
Particles are thereby produced through the inflaton's coupling to other matter fields. If the 
matter fields are dynamical and chaos reigns, then the process of entropy production would 
deserve rethinking. The importance of parametric resonances, which are closely related to 
chaotic behaviour, has already been stressed in this context |20J21[] . 



It seems appropriate to consider how chaotic dynamics might impact on "chaotic" infla- 



tion ||17|| . In chaotic inflation different patches of the universe are taken to have different 
values of the inflaton and matter fields. In some patch, it is argued, the inflaton is suffi- 
ciently high up in the potential and the matter fields are sufficiently small so as to permit 
a long-lived inflationary epoch 0. For initial field values deep within a basin, away from 
the fractal borders, the usual arguments hold and the chaotic inflation paradigm is largely 
unaffected. However, if in a given patch the field values are near a fractal basin boundary, 
it can become difficult to find a patch of any size across which the conditions are regular 
enough to allow this thinking. Even the slightest variation in the initial conditions across 
the patch will lead to an entirely different outcome. Due to the self-similar nature of the 
fractal, no matter how small you try to make the patch, there will still be slight variations 
in the conditions and hence the outcome. For cosmological conditions in the vicinity of the 
fractal basin boundary then, the simple FRW thinking must be abandoned. 

For similar reasons, caution would be needed for slow-roll initial conditions as well. 
In fact the slow-roll scenario will likely be more fragile as the inflaton is more easily kicked 
around. For chaotic initial conditions by contrast the field high up in its potential is resiliant 
against the influence of kicks and bumps. 

Aside from direct observational effects there are also some important theoretical im- 
plications raised by chaotic evolution. For example, chaotic systems are characterised by 
an entropy, the Kolmogorov- Sinai entropy, which is related to the spectrum of Lyapunov 
exponents. This introduces a chaotic arrow of time in addition to the cosmological and 
thermodynamic arrows of time. As well, it raises the question of a possible connection be- 
tween the Kolmogorov- Sinai entropy and the thermodynamic entropy released at the end of 
inflation. Another issue raised by chaotic dynamics concerns the recovery of a semiclassical 
limit in quantum cosmology due to the breakdown of the WKB approximation in chaotic 
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systems [^2 . 

We have suggested a few implications for chaotic dynamics in largely unexplored terrain. 
Chaos theory grew out of Poincare's study of the solar system, and with its development 
came insights into the intricate structures of our neighborhood such as the asteroid belt and 
Saturn's rings. It seems fitting for chaos to have an impact not only on the evolution of the 
solar system but also on the birth of the universe. We are left to ask if chaotic fingerprints 
have been left on the large-scale landscape as they were on the landscape of our own solar 
system. 
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FIG. 1. Phase space trajectories in the (a, a') and (ip,ip') planes for universes with 
non- interacting scalar fields. The solid line is the separatrix and the cross hatched regions mark 
the big crunch basin of attraction. 
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FIG. 2. The correlation between initial and final values of the scalar field (a) is compared to 
the apparently chaotic behaviour seen in (b) where the sampling rate is ten times lower. 



FIG. 3. The three possible outcomes for the universe. In each case the solid line is the scale 
factor a and the dashed line is the scaled acceleration da 3 . The initial values for are —23.31, 
—23.32 and —23.33 respectively. 
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V>1 

FIG. 4. The basins of attraction for universes similar to those shown in Fig. 3. 



FIG. 5. A detail of Fig. 4. 
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FIG. 6. Finding the fractal dimension for Fig. 5. The solid line is a least-squares fit to the 
box counting data. The dimension was found to be ds = 1-58 ± 0.02. 



a' 



FIG. 7. A slice in the a - a' plane which intersects the ipi - ip[ plane of Fig. 4. along the line 

y>i = l.o. 
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V>1 

FIG. 8. The road to chaos: As A is incremented from 0.5 to 2.0, the nonlinear distortion of the 
attractor basin boundaries mounts. Once A exceeds 1.0 the mixing is so strong that the boundaries 
become fractured, and eventually, fractal. The graphs were generated for the choice of parameters 
and initial conditions A = 0.00002, a = 10, ^2 = 5.0, ij/ 2 = 10 -°i m i = °- 2 > m 2 = 0.1. 



FIG. 9. The correlation between ipu (vertical axis) and ipi m (horizontal axis) on a ip[ = slice 
through the big crunch basin of attraction of Fig. 8. (A = 2.0). 
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FIG. 10. Basins of attraction in the <J>i - $i plane for universes containing two minimally 
coupled scalar fields. 



$1 

FIG. 11. A detail of Fig. 10. where the dimension is ds = 1-54 ± 0.02. 
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FIG. 12. (a) A trajectory spiraling into the de Sitter attractor. (b) A nearby trajectory which 
flows out to the big crunch attractor. 



FIG. 13. The basins of attraction in the a - a plane for universes containing a minimally 
coupled scalar field and a massless, conformally coupled scalar field. 
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FIG. 14. A detail of Fig. 13. where the dimension is d B = 1.484 ± 0.005. 
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